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■ Using certain finite-dimensional stable range of the nonlinear terms, we obtain large families 
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of exact solutions parameterized by functions for the equation of nonstationary transonic gas 
flows discovered by Lin, Reisner and Tsien, and its three-dimensional generalization. 



a^ : 1 Introduction 

00 

! Lin, Reisner and Tsien [LRT] found the equation 

o , 

^ . '^Utx + U^U^x - Uyy = 0. (1.1) 

- ^ , for two-dimensional non-steady motion of a slender body in a compressible fluid, which 
' was later called the "equation of nonstationary transonic gas flows" (cf. [Ml]). Mamontov 
obtained the Lie point symmetries of the above equation in [Ml] and solved the problem 
of existence of analytic solutions in [M2] . The three-dimensional generalization: 

2Utx + UxUxX - Uyy - U;,Z = (l-2) 

was studied by Kucharczyk [K] and by Sukhinin [Ss]. Indeed, the Lie point symmetries 
of the equation (1.2) were found in their works. Sevost'janov [Sg] found explicit solutions 
of the equation (1.1), describing nonstationary transonic flows in plane nozzles. 
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In this paper, we present a new approach based on the fact that the nonhnear terms 
keep some finite-dimensional polynomial space in x stable. We obtain a family of solutions 
of the equation (1.1) blowing up on a moving line y — f{t), which reflect partial phenom- 
ena of gust, and a family of smooth solutions parameterized by six smooth functions of 
t. Moreover, we flnd a family of solutions of the equation (1.2) blowing up on a rotating 
and translating plane cosQ;(t) y + sinQ;(t) z = fit), which reflect partial phenomena of 
turbulence, and a family of solutions polynomial in x parameterized by time- dependent 
harmonic functions in y and z, whose special cases are smooth solutions. In particular, 
we find all the solutions polynomial in x and y for the equation (1.1) and all the solutions 
polynomial in x. y and z for the equation (1.2). Since our solutions contain parameter 
functions, it can be used to solve certain boundary-value problems for these equations. 

Lie group method is one of most important ways of solving differential equations. 
However, the method only enables one to obtain certain special solutions. It is desirable 
to find more effective ways of solving differential equations. Indeed, we do find one of 
solving the above nonlinear partial differential equations. 

On Mamonotov's hst of the Lie point symmetries of the equation (1.1) (e.g. cf. Page 
296 in [I]), the most sophisticated ones are those with respect to the following vector 
fields: 



Xi = 3a{t)dt + {a'{t)x + a' {t)y^)d^ + 2a'{t)yd, 

-a'{t)u + a"{t)x^ + 2a"\t)xy^ + -a^^\t)y'' 

3 



y 



+ 



(1.3) 



X2^l3'{t)yd, + I3{t)dy + 



2(5"{t)xy+'^-P"\t)y^ 



du, (1.4) 



^3 = + [27'(^)^ + 2YW]du, (1.5) 

where a, (5 and 7 are arbitrary functions of t. Among the known Lie point symmetries 
of (1.2) (e.g. cf. Page 298 in [I]), the most interesting ones are those with respect to the 
following vector fields: 

Xi = y^'^dt + {^x + ^(y^ + -2^)^ 9x + ^tydy + Ztzd, + [x^ - 3tu)du, (1.6) 

X5 = g\t)yd, + 9{t)dy + [2g"{t)xy + g"\t) (^^ + yz^^ ^ d^, (1.7) 

X, = h\t)zd, + hd, + [2h"{t)xz + g"\t) (^y + zy^^ ^ d^, (1.8) 
X, = ad, + [2<T'{t)x + (T"{t){y'' + z')9„, (1.9) 
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where g, h, a and a are arbitrary functions of t. Ryzhov and Shefter [RS] found some 
invariant solutions of the equation (1.2), which represent time-dependent flows in a circular 
or plane Laval nozzle. 

First we find that the group invariant solutions with respect to the vector fields Xi- 
Xr are polynomial in x with degree < 3. Then we examine the equations (1.1), (1.2) 
more closely and observe that this phenomena is essentially caused by the fact that the 
nonlinear term u^u^x keep the following polynomial subspace stable: 

R + Mx + Rx^ + Rx^, (1.10) 

where R stands for the field of real numbers. This observation suggests us a new ansatze 
of solving the equations (1.1) and (1.2). Since the equation (1.2) contains the Laplace 
operator 9^ + 9^, our approach to (1.2) will involve harmonic analysis and sophisticated 
integrations. For simplicity, we will solve the equation (1.1) in Section 1.2 although it can 
be viewed as a special case of (1.2). Exact solutions of (1.2) will be found in Section 3. 

2 Two-Dimensional Case 

In this section, we study solutions polynomial in x for the Lin-Reisner-Tsien equation 
(1.1). By comparing the terms of highest degree in x, we find that such a solution must 
be of the form: 

n = fit, y) + g{t, y)x + h{t, y)x' + i{t, y)x\ (2.1) 

where f{t,y), g{t,y), h{t,y) and ^(t,y) are suitably-differentiable functions to be deter- 
mined. Note 

Ux^g + 2hx + 3^x^ ^2h + 6^x, (2.2) 

'^tX = + 2htX + 3^tX^, Uyy = fyy -\- g yyX -\- kyyX^ + ^yyX^, (2.3) 

Now (1.1) becomes 

2{gt + 2htx + 3^tx^) + {g + 2hx + 3^x^) {2h + 6^x) - fyy - gyyX - hyyX^ - iyyx^ = 0, (2.4) 
which is equivalent to the following systems of partial differential equations: 

iyy = 18e^ (2.5) 

V = 6^^ + 18^/1, (2.6) 
= Aht + Ah^ + Qig, (2.7) 
fyy^2gt + 2gh. (2.8) 



hy = ^lw ^ /2.10) 



First we observe that 

C = r- ^ (2-9) 

is a solution of the equation (2.5) for any differentiable function j3 oit. Substituting (2.9) 
into (2.6), we get 

12p'{t) , 18 
Denote by Z the ring of integers. Write 

h{t, y)^Yl «^(^)(^?/ + Pit))'- (2-11) 

Then 

hyy = E + 2)(^ + l)ai+2mVSy + mf- (2-12) 
Substituting (2.11) and (2.12) into (2.10), we have 

g 3[(i + 2)(z + 1) - 6]a,+2mV3y + mY = " (2-13) 

equivalently, 

-12a_i(t) = -12/3'(t), 3(i + 4)(i-l)ai+2(t) =0, i 7^ -3. (2.14) 

Thus 

where a and 7 are arbitrary differentiable functions of t. 
Note 

^ ^ 2«(t)/3^(0 ^ a'(i) /3'(i)2 



iV3y + m)' iV^y + m)' iVSy + m)' 

+ + 37(t)/3'(t)(V3|/ + /?(t))2 + 7'(t)(x/3y + /^(t))^ (2.16) 



and 



(V3|/ + /3(t))4 (v/3t/ + /3(t))3 (v/3y + /?(t))2 
+2a(i)7(i)(V3y + m) + 2j{t)P'{t){V3y + m)' + l{t)'{^y + m)'- (2-17) 

Substituting the above two equations into (2.7), we have: 

6 4a' Ap" 



+8aj{V3y + /?) + 20-fp'{V3y + (5f + Ai{V^y + (5f + A-t^{V^y + /?)^ (2.18) 



Write 

g{t,y)^Y.^i{t){V3y + py. (2.19) 

Then 

9yy = E + 2)(^ + l)k+2{t){V3y + (3)'. (2.20) 
Substituting (2.19) and (2.20) into (2.18), we get 

3(i + 3)z6i+2 - 0, ^ 7^ -4, -2, -1,1, 2, 3, 6, (2.21) 

h 2/?" 2«7 

0-2 = y , Oo = — ^, Oi = — ^, h = (2.22) 

2/3'7 , 2y 272 
&4 = ^, h^-, bs^—. (2.23) 

Therefore 

+ - ^ - ^(VSy + /3) + p(V3|/ + /?)2 
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3{VSy + (3f VSy + (3 3 3 



+^{V3y + + ^(V3|/ + + 1^(731/ + + ^(V3|/ + /?)«, (2.24) 

where a and p are arbitrary differentiable functions of t. 
Observe that 

^ _ 2a^p' ^ {2aa'-3af3') ^ a' 2a" 2(5' (5" 



3(v^y + /3)3 3(V3y + /3)2 V3y + /3 3 
^2(3p^^-^^^^ + /3) + (p' + 2^7/30(^31/ + 



^ 2.'7 + 2.y + 8(/30^ ^^^ ^ ^ 18/3-7 H- 28/3Y ^^^ ^ 
+^(V3|/ + /3)^ + ^(V3y + (5y + ^^{^2,y + (5)\ (2.25) 



27 V ^ 81 ' ^ 81 

3^^ ^ ^-^2^/ 3^/^ _ 2ctct' _ 2 (a/3" + a'/3') 
^ " 3( V3|/ + /3)4 ^ 3( VSy + /3)3 ^ 3(^y + /3)2 ~ 3( VSy + /3) 

+ap - ^ + (.^ + /3'p)(V3, + /3) + + 
^ 18(/30-7-16.7- ^^ + /3)3 + ^^^^^^(73, + ,Y + 7P(/% + /3)^ 

+^(^. + /3)« + ^(^^ + /^)^ + !r(^^ + (2-2^) 
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Substituting (2.25) and (2.26) into (2.8), we obtain 

_ 2a^ Qaa - 2a^(3' 6a' - A{aj3" + a' 13') 

~ 3(x/3y + /3)4 ^ 3(V3|/ + /3)3 ^ 3(A/3y + /3) 

+2.p-if::-gf::+ ^-^^+^f^-^^" (V3.+/3) 

^20./3-7 + 67. + 6p'^^^ ^ ^ 18O(/30^7 + W + 36.-7 ^ 
3 27 



+^(V3, + + ^(V^!/ + «^ + 1^ (V3y + (2.27) 



Thus 



- 4(a/3" + a'/3')](A/3y + P)MV3y + P)-1] + ^^^^^^^§f^^{V3y + Pf 

+ 



54 V ^ ' / 4Q5 



^1 ' ^ ' ' 1701 ^ ^ ' ^ 243 

+ + .^{^y + (5Y^ + ^(V3y + (2.28) 
243 ^ y 10935^ ^ ^ 9477^ y a-v > v ; 

where ^ and '& are arbitrary functions of t. 

Theorem 2.1. W^e have the following solution of the equation (1.1) blowing up on 
the surface \/?>y + P{t) — 0; 

u = — ^ \ ^ h ^= h 7(V 3|/ + /?) V 

«^ 2q;' 2(3'' r m / /- m2 

+ 27(v5!, + /3P+5(7S^ + ' + *^"''^ 1 " 

+h(M,' - 4(«/3" + a'/3')](V3!, + /J)[ln(v^!, + /J) - 1] + 3a'7 + 9/?> - 2,3-' ^ 

+^^(^^ + m9 + + + ^(x/Sy + (2.29) 

243 ^ 10935^ 9477^ ' ^ ' 

where a, /3, 7, cr, p, ^ and '& are arbitrary functions oft, whose derivatives appeared in the 
above exist in a certain open set ofR. 

When a = j = a = p = 9 = 'd = 0, the above solution becomes 

= + —iVSy + (3)x - —^y^ - —^iV^y + Pf. 2.30 

(V3y + /3)2 v^y + /3 3 ^ ^ 3 ^ 27 ^ ^ ^ ' 



Take the trivial solution ^ = of (2.5), which is the only solution polynomial in y. 
Then (2.6) and (2.7) become 

hyy = 0, gyy = Aht + (2.31) 

Thus 

h = a{t) + i3{t)y. (2.32) 

Hence 

Qyy = 4q;^ + Aa' + 4(/3' + 2a(3)y + 4(3^. (2.33) 

So 

^ = ^ + at/ + 2ia' + a')y' + ^(/9' + 2«/5)|/3 + i/?V, (2.34) 
where 7 and a are arbitrary functions of t. Now (2.8) yields 

fyy = 2{a-)- + -'/')+ 2{aa + + a')y + 2{pa + 2a^ + 6aa' + 2a")y^ 

20^2^ + 12a/5' + 20a'/? + 4/3" , 10aP^ + 8(3(3' 4 2 ^3 , 
+ ^ 3 Z+g/gV- (2.35) 

Therefore, 

/ = r + py + {a^ + -f')y^ + ^ {aa + /?7 + cr')?/^ + ^ (/?(7 + 2a^ + 6aa' + 2a")y^ 

, 5a'P + 3aP' + 5a'P + P" , ha(3^ + A(3(3' e ^ 1 ^3 7 to 

H y H y ^ (> y ■ 2.36 

15 ^ 45 ^ 63 ^ ^ 



Theorem 2.2. The following is a solution of the equation (1.1): 
u — (a + (3y)x^ + 



7 + ay + 2(^2 + a')y' + \{(3' + 2a/3)y3 + 



X 



+r + py + (q;7 + 7')^^ + ^(acr + /37 + (7')y^ + ^(/?(7 + 2q;^ + 6aa' + 2a")y^ 

o D 







63 



(2.37) 



15 " 45 

where a, P, 7, a, p and r are arbitrary functions of t, whose derivatives appeared in the 
above exist in a certain open set of R. Moreover, any solution polynomial in x and y of 
(1-1) must be of the above form. The above solution is smooth ( analytic) if all a, (3, ■j, a, p 
and T are smooth (analytic) functions oft. 

Remcirk 2.3. In addition to the nonzero solution (2.9) of the equation (2.5), the other 
nonzero solutions are of the form 



e = p.(V3y + /3(i)), 



(2.38) 



where is the Weierstrass's elhptic function such that 

p'Xwf = A{p,{wf-i), (2.39) 

and t is a nonzero constant and (5 is any function of t. When (5 is not a constant, the 
solutions of (2.6)-(2.8) are extremely comphcated. If (3 is constant, we can take /3 = by 
adjusting u. Any solution of (2.6)-(2.8) with /i 7^ is also very comphcated. Thus the 
only simple solution of the equation (1.1) in this case is 

u = p,(v^y) x\ (2.40) 

3 Three-Dimensional Case 

By comparing the terms of highest degree, we find that a solution polynomial in x of the 
equation (1.2) must be of the form: 

u = f{t, y, z) + g{t, y, z)x + h{t, y, z)x^ + ^(t, t/, z)x^, (3.1) 

where f(t,y,z), g{t,y,z), h{t,y,z) and ^{t,y,z) are suitably-differentiable functions to 
be determined. Note 

Ux — g + + 3^x^, Uxx — 2h + 6^x, (3-2) 

"l^tX — fl't + 2/l(X + 3^(3;^, Uyy — fyy + g'^/J/^ "I" ^J/J/"^^ ^2/2/"^^' (^'S) 

Uzz = fzz + QzzX + KzX^ + ^zzX^. (3.4) 

Now (1.2) becomes 

2{gt + 2htx + 36a:') + {g + 2hx + 3^x^){2h + 6^x) - {fyy + U) 

-i9yy + 9zz)x - {hyy + h^z)x'^ - {^yy + ^^^)a;^ = 0, (3.5) 

which is equivalent to the following systems of partial differential equations: 

Cyy + Cz.^lSe, (3.6) 

hyy + h,, = 6^t + IS^/i, (3.7) 
9yy + 9zz = + + 6Cg, (3.8) 
fyy + = 2hgt + 2gh. (3.9) 

First we observe that 

1 



( V3(y cos a{t) + z sin Q;(i)) + P{t)y 
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(3.10) 



is a solution of the equation (3.6), where a and (3 are suitable differentiable functions of 
t. With the above ^, (3.7) becomes 

12(\/3q;'(— l/sincK + zcosa) + /9') 18/i (ot^\ 
(v3(?/cosq; + 2; sin a) + /?)-^ (v 3(l/ cos a + 2; sin a) + 

In order to solve (3.11), we change variables: 

C = •\/3(cosa 7/ + sina 2;) + /3, 77 = \/3(— sina 7/ + cosct z). (3-12) 

Then 

COS adi; — sin a dj,) , dz — Vs(sm ad,^ + cos a dj,) . (3.13) 

Thus 

dy + dl = 3(cosq; — since 5^)^ + 3(sinQ; + cos a 9^)^ = 3(9^ + 9^). (3-14) 

Note 

dtiO^a'rj + P', dt{rj)^a'{P-C). (3.15) 
The equation (3.11) can be rewritten as: 

kc + Kn = -4(a'^? + /5')C"^ + 6C"'/i- (3.16) 
In order to solve the above equation, we assume 

h = J2<^iit^ri)C- (3.17) 

Now (3.16) becomes 

Yl ^imC + E - = + ^')C-' + 6 a,C-^ (3.18) 
which is equivalent to 

a_3^^ + 2a_i = -4(a'r7 + /3') + 6a_i, a^^,, + (i + 2)(i + l)ai+2 = 60^+2 (3.19) 
for — 3 7^ i e Z. Hence 

a_i = ^a-3r?r) + Q;'?7 + /?', (i + 4)(i - 1)0^+2 = -a^^r? for - 3 ^ i e Z. (3.20) 

When i = — 4 and i = 1, we get a_4^^ = ai^^ = 0. Moreover, a_2 and 03 can be any 
functions. 
Take 

oa = (T, a_2 = p, a_i = a'r? + (3-21) 



Oi = a-i-2i = a-2-2i = for < i e Z (3.22) 

in order to avoid infinite number of negative powers of C in (3.17), where a and p are are 
functions of t and 77 differentiable in a certain domain. By (3.20), 

~ nS^^T^ " (2fc + 5)(2fc + 3)(2^ + l)!' ^^-^^^ 

^-^^ ^ nl.(20(2^-5) = w ■ ^'-'^^ 

Therefore, 

h . (.w/>')c- + E(-i) ^,....jiiu. 



fc=0 



■(2A; + 5)(2A; + 3)(2A; + 1)! 



{2k-l){2k-3)df{p)C-' 
+ 3(2^)! ^^-^^^ 

is a solution of (3.16). 

By (3.12) and (3.14), (3.8) is equivalent to 

^CC + 9m = + ^h^ + 2C^g. (3.26) 

Note 

00 

h, = {a"v + P" + {a'fPX-' - {a'f - {a'v + PfC' + ^(-l)'^''^ 

k=0 



X 



(2A; + 5)(2A; + 3)(2A; + 1)! 3(2A;)! 



. , .nr IH'W , (2fe-l)(2fe-2)(2fe-3)af(p)C-^ , 

+ ^ + ^ ^[ (2fe + 5)(2fe+l)! + 3M! ^ ^ ^ 

For convenience of solving the equation (3.26), we denote 

A A °° 

^ht + -h'=J2Ht,v)C (3.28) 

i=-4 

by (3.25) and (3.27). In particular, 

b-4 = ^p', &-3 = 0, (3.29) 

4 4 

&-2 = g(pt + a'f^Pv) + gPrmP, (3.30) 

6_i = ^K'r^ + + {a'fp - a'pr,] + ^{a'r^ + /3')p,^, (3.31) 

^-0 = + g (Pt.. + «'/3p,,,) + -aj(p)p + -pj^. (3.32) 
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Suppose that 

g = Y,c,{t,r^)C (3.33) 

is a solution (3.26). Then 

oo 

ieZ r-=-4 ieZ 

equivalently 

(i + 3)iQ+2 = fcj - (^ + 3)rc^+2 = -Crr^n, r < -A<i. (3.35) 

By the second equation in (3.35), we take 

Cr^O for r < -4 (3.36) 

to avoid infinite number of negative powers of C in (3.33). Letting i = —3, 0, we get 

6_3 = c_3^^, bo — Cqjjjj. (3.37) 

The first equation is naturally satisfied because c_3 = — c_5^^/10 = 0. Taking i — 
—2, —4 and r = — 6 in (3.35), we obtain 

Co = ^C-2r,r, " ^&-2, C_2 = ^6-4. (3.38) 



So 



Co = Id^ib.,) - h^. (3.39) 



Thus we get a constraint: 

bo = Idiib-d - ^9,^(6-2), (3.40) 



equivalently. 



^ + g(PtW + a f3pr,nri) + g<9^(p)p + T^Pr^n 



3 9'^ " 9"'?^'^'^ 27 

i4/„2 



Thus 

48(p,,, + a'/3p,,,) + 69;;(p)p + 2pJ^ - 9d'^{p') + I2d'^{p,,p) = 72(a')'. (3.42) 

It can be proved by considering the terms of highest degree that any solution of (3.42) 
polynomial in rj must be of the form 

P = 7o(^)+ 71(^)^ + 72(^)^7'. (3.43) 
11 



Then (3.42) becomes 



127^ - 2072^ = 9(a') ■ 



(3.44) 



So 



2e 



OL = —\l 372 — 57I ^ a 



2e 



37^ - 57|cit, (3.45) 



where e = ±1. Replace /3 by —j3 if necessary, we can take e = 1. Under the assumption 
(3.43), 



fe=0 



(2A; + 5)(2A; + 3)(2A; + 1)! 



and 



6_i = ^[a"rj + P" + {a'fp - a'p,] + ^{a'rj + /3')72, 
4(0;')' 4 , 4 2 



(3.46) 

(3.47) 
(3.48) 
(3.49) 



by (3.30)-(3.32). 
Denote 



For any real function F{t, rf) analytic at 77 = 770, we define 



r=0 



Note 



15af(a)C 



k=0 



(2A; + 5)(2A; + 3)(2A; + 1)! 



15C' 



•^0 Vfc=o 



(a)rf 



(2A; + 5)(2A; + 3)(2A;)! 



(3.50) 



(3.51) 



15 / r2 / 

/o ^0 



A:=0 

= 15C-^ / r3 / rj 
'0 -/o -/o 



(2A; + 5)(2A;)! 



.A;=0 



dri dT2 drs 



= — C"^ / '^3 / '^2 / [cr(t, ^7 + \/^Ti) + a{t, 7] - ^/^Tl)]dTl dT2 dr- 
^ Jo Jo Jo 

00 

E(-i)'= 

.fe=0 



(3.52) 



d. 



159f (fT)C3+2fc 
(2/e + 5)(2/e + 3)(2A; + 1)! 



15af+i(a)C'+2*^ 



(2A; + 5)(2A; + 3)(2A; + 1)! ' {2k + 5){2k + 3){2k + 1)\ 

^{2k + b){2k + ?,){2k + iy.^ '^^^^ ^ 



k=0 



k=0 



(2A; + 5)(2A;+ 1)! 
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= -TtC^ I ^3/ T2 I [(7t(i,?7 + V^Ti) +(7t(i,77- V^Ti)]dTidT2dT3 

^ JO -/o -/o 

+ 0/-2 / ^2 / Ti[a(i, 77 + ^Ti) - (7(i, 77 - V^Ti)]dTi 

Jo Jo 



+^{a'r) + /3')C"' / T-| / 77 + V^n) + (7(i, 77 - y/^ri)]dn dr2. (3.53) 
^ Jo Jo 

According to (3.28) and (3.50), we have 



+10C~^ / '^3 / T2 / [o-t(^, yi + \/^7-i) + (Tt(t, T] - v^ri)](iri rfrs rfrg 
Jo Jo Jo 

+ ^^-^^ / / ri[a(t, 77 + \/^ri) - (7(t, y] - x/^ri)]o?ri 

s Jo Jo 

+ 10(a'7/ + /?')C^' I r [^(^' ^ + V^n) + ^(^' ^ - V^Ti)]dTi dT2 

Jo Jo 

J T2 j [(7{t, 77 + ^^-^Tx) + (7(i, 77 - V^Ti)]dTi ^Tg j 

xC-^(pr'+(a'77 + /3')r' + |). (3.54) 
Now 

C-2 = y (3.55) 

by (3.29) and (3.38). According to (3.35) with i — —3, 0, c_i and C2 can be arbitrary. For 
convenience, we redenote 

C-i = K{t,r]), C2 = u}{t,r]). (3.56) 

Moreover, (3.29), (3.39) and (3.47) imply 

2 2 
Co = y - g(Pt + a'/^Pr,) + g72P- (3.57) 

Furthermore, (3.31) and (3.35), 

ci = ^ - ^K'77 + P" + {a'f(5 - a'p,\ - ^(a'77 + /3')72. (3.58) 

In addition, (3.35) and (3.56) yield 

{-ir^df-\.) f (-i)^-(.+i)(2.)! 

'^'=+^-2()t + 2)(2)t + 2)!+^ ()t + 2)(2)t + 2)! ^^"^^^ 

_ (-l)^.i3af+2(,) f (-l)^-(2. + 3)(2. + !)! (,_, 

- (2A; + 5)(2A; + 3)! + ^ (2A; + 5)(2A; + 3)! ^^''^'^ ^^'^^^ 



for < A; e Z. 
Set 



^ oo 



i=3 



,fc=0 



(2fc)! 



3af(a;)C2'= 



k=0 



{2k + 3){2k + 1)\ 



oo 



fc=0 i=0 

oo k 

k=Q 1=0 



+ l)(2i)! 
(fc + 2)(2A; + 2)! 



2k+3 



l)'^-*(2i + 3)(2i + l)! 
(2A; + 5)(2A; + 3)! 



2fc+4 



Note 



39f(c.)C^'= 



^ > i2k + 3){2k + i)\ 

-C"^ / T2 / [a;(i, 77 + v^Ti) + rj - V^Ti)]dTi dT2. 
^ Jo Jo 



Moreover, 



XT; m n 9^(*(^,p,.>)(t,77,0) 



for < i e Z. 



Thus 



^{P,p,a,K,iv) (t, V, C) = 7;C~^ / ^2 / + V^n) + uj{t, rj - ^/^n 

^ Jo Jo 

h>h {2i + l){k + 2){2k + 2)\ ^ 

^ (-l)^-(2. + 3)9f-^c}f+^(v&(,.,..))(t,r7,0) 
hh {2t + 2){2k + 5){2k + 3)\ ^ ' 

in which the summations are finite if a{t, rj) is polynomial in rj. Now 



6 2 



Jo Jo Jo 



X 



c? (? 2 2 
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Denote 

A(/3,p,a,K,a;>(i, V, C) = ^(^^^ + dt{g)). (3.67) 
Then the equation (3.9) becomes 

fvv + fa = ^{P,P,^,K,u,) {t,V,0- (3.68) 



Set 



Then 



Hence 



w^^-±^, W^^-^^. (3.69) 



drj = \{d^ + d^), dc = ^{d^-d^). (3.70) 



d^ + dl^d^d^. (3.71) 

A complex function 



G{iJ,) is called bar-homomorphic ii G{ijl) — G{ijl). (3.72) 

For instance, trigonometric functions, polynomials with real coefficients and elliptic func- 
tions with bar-invariant periods are bar-homomorphic functions. The extended function 
F{t,n) in (3.51) is bar-homomorphic in /i. Now we have 

f ^ /_A(/3,p,<x,«:,c^>(i,/^ + 7^, a/^(7I-/^))c?7^(^/^ + x(^,?7 + a/-1C) 

•J wo •JWq 

+x{t, V - V^C) + y/^Ht, V + V^C) - ^{t, V - V^C)], (3.73) 



where wq is a fixed complex number and x(t,yu), i'{t,fi) are complex functions in real 
variable t and bar-homomorphic in complex variable /j,. 

Theorem 3.1. In terms of the notions in (3.12), the function u = ("'^x^ + hx'^+gx + f 
is a solution of the equation (1.2) blowing up on the hypersurface ^/3{cos ay +sm a z)+ (3 = 
0, with h given in (3.65), g given in (3.66) via (3.54) '^''^d (3.64), ^^f^d f given in (3.13) 
via (3.67). The involved parametric functions p is given in (3.43), a is given in (3.45), 
and (T, uj are real functions in real variable t and rj. 

Next we want to find a more explicit formula when (T = fi; = a; = x = i/ = 0. In this 
case, 

h = pC' + {a'rj + (5')C' + I (3.74) 



and 

9 = ^C' + §-l{Pt + a'Pp,) + lj,p 

A[a"rj + (3" + {a'f (3 - a'p,]C - \{a'r) + /3')72C- (3-75) 
Moreover, (3.15) and (3.43) yield 

-6p, - Qa'Pp, - 272P) - \p{a"r^ + + («')'/5)]C-^ ^ ^ 2|p 

-giPtt + 2q; /^(ptr, + (a P) Pr, + 2(q; (3) 72) + + h + — 

-^-{a'r^ + /30(a"r; + P" + (a')'/9 - «'p.) - ^72(a'r; + /3')^ + 
+ |(19(aO^/3 - 5a"r, - 5/3") - - ^^^l{a'r, + /3') - ^(a'"r; + P'" 

+ {arP' - a"p, + 3aV'/3) - ^(a'7^,7 + /J^DIC + - ^^^]C'- (3-76) 

By (3.67) and (3.68), we can find /. 

Corollary 3.2. In terms of the notions in (3.12), we have the following solution of 
the equation (1-2) which blows up on the hypersurface \/3(cosa y + sin a z) + 13 = 0; 



u = ("V + 



PC' + {a'v + P')C-' + ^ + i^c' + ^ - 9(p, + a' Pp.) 



6 



3 ' 3 3 



p^ _ {o^rj+fy ^121 72p|+2^ 

27^ 9 ^ 27 ^ 54 ^ ^ ^ ^ 

+ ^C'(12 InC - 25) + i^a'pp, + ^^^^|^(9p; - 6p, - Qa'Pp, + A^,p) 
~p{a"rj + P" + {a'fpminC - 1) - ^[29a'72P, + {a'v + P'){20^i - H) 
-%a'ff3l2 - 3{a'p,, + ^,{a"r, + (3") + a'\)]C^(6 InC - 11) + + ^ 

>a ^( >a\> >a\2 \ , 72Pt , 13a'/372P^ 72P^ 
- g(p« + 2q; (3ptr^ + (a /3) + 2{a (3) 72) + + + 

+^ - ^(«'^ + /90(«^ + + - «'P.) - ^72(a'r? + /?')'](' 

^y-^2rr 81 ^ ) /3 - 5a - 5/5 ) 243^2(" V + P ) 

-^(o^"'v + P'" + (o^fP' - o^"Pv + 3«V'/3) - l(a'7^77 + /^^DlC' 

r7272 72' , 77I , 4(aOS i^4 /o 

^"6" 27 + 486 ~9~ + ^^^^' ^'^"^^^ 
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where p is a third-order dijferentiable function as in (3.4.3), a is given in (3.45) and (5 is 
any third-order dijferentiable function of t. 

Finally, ^ = is the only solution of the equation (3.6) polynomial in y and z. Under 
this assumption, the equations (3.7) and (3.8) becomes 

hyy + Kz = 0, Qyy + Qzz = + 4/i^ (3.78) 

The first equation is a Laplace equation whose solutions are called harmonic functions. 
It can be proved as (3.25) by power series that the general solution of the first equation 
is: 

h^{a + y/^p) it,y + y/^z) + (a - ^/^p) {t,y - y/^z) , (3.79) 

where (j(t,fi) and p{t,p) are complex functions in real variable t and bar-homomorphic 
in complex variable p (cf. (3.72)). Set 

w — y -\- V^z, w — y — y/^z. (3.80) 

Then the Laplace operator 

d^y + dl^Ad^d^. (3.81) 
The second equation in (3.78) is equivalent to: 

dnjdwig) ^ht + h'^ ^ {at + V^pt) (t, w) + {at - y/^pt) {t, w) 

+ [{a + V^p){t,w) + {a- V^p){t,w)]^. (2.82) 

Hence the general solution of the second equation in (3.78) is: 

9=1 {{'^t + V^Pt){t:Pl) + {at - V^pt){t:TH) 

J w\ J wi 

+[{a + V^p){t,pi) + {a- V^p){t,-iM)]'^}dpi djM 

+{K + V^u;){t,w) + {K-V^u;){t,w), (3.83) 

where K{t,p) and uj{t,p) are complex functions in real variable t and bar-homomorphic 
in complex variable p (cf. (3.72)), and wi is a complex constant. Furthermore, (3.9) 
becomes 

dMf) = li9t + 9h). (3.84) 

Observe 

pw pw 

9t= _ [2{{a + V^p){t,p,) + {a-^/^p){t,pl)){{at + ^/^Pt){t,Pl) 

+ (crt - ^^pt)(t,7Il)) + {att + \/^ptt){t,Pi) + {att - \f^Ptt){t,W)]dpi dpi 

+(Kt + %/^Ut){t,w) + {Kt - V^u;t){t,w). (3.85) 
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Thus the general solution of the equation (3.9) is: 

x(((jt + ^/^pt){t,iii) + (cTj - ^/^pt){t,-fl{)) + {utt + ^/^ptt){t,^^l) 

+ (crtt - v/^ptt)(t,7ir)](i/xi 6411+ {Kt + v^cc;t)(t,/U2) + («t - ^^a;t)(t,7I^) 

+ [((7+ ^p)(t,/X2) + (a - v^p)(i,7I^)][(/^+ V^C^)(t,/X2) 

+{K-^/^uj){t,iI^)+ I / {((7t + v^pO(^,Afi) + (<^t- v^pO(^,Afr) 

J UIl J Ul 

+ [((7 + \/^p)(t,/ii) + (a - v^p)(t,7^)]^}(i/xi (i7Zr]}(i/X2 ci/z^ 

+(X + w) + (x - w), (3.86) 

where xi't^p) and i^{t,p) are complex functions in real variable t and bar-homomorphic 
in complex variable and W2 is a complex constant. 

Theorem 3.3. In terms of the notions in (3.79), the following is a solution polynomial 
in X of the equation (1.2): 

u=[ia + ^/^p){t,w) + {a-^/^p){t,w)]x^ + { I / {(cx^ + v^p)(t, /xi) 

J Wl J Wl 

+{at - V^p){t,jH) + [{a + V^p){t,iii) + {a- v^p)(t,7IT)]2}d/ii djH 

+ {K+V^Lj){t,w) + {K-V^u){t,w)}x+- I { _ / 

J W2 J W2 J Wl J Wl 

[2{{a + V^p){t,i^i) + {a- V^p){t,jHMat + V^pt){t,i^i) 

+{at - y/^pt)it, /TT)) + {(Ttt + y/-lptt)it, l^i) + {ou - \f-ipu){t, A*T)](^/ii rf/TI 

+ (fi;t + \/^ujt){t, p,2) + {Kt - \/~^uJt){t,Jpi) + [{a + yf-[p){t, 1x2) 
+ ((T - ^^^p) {t, -p^)] [{k + ^/^uJ) {t, 1^2) + V^ou) {t, -p^) 

+ /_ / {i<^t + ^Pt)it, 1^1) + {(Jt - ^Pt)it, in) 

J Wl J Wl 

+ [((T + v^p)(t,/ii) + (0- - ^/^p){t,JIl)f}dni djH]}dii2d-p^ 

+{X + V^'^){t:w) + (x - V^u){t,w), (3.87) 

where a{t, p), p(t,p,), n(t,p,), u!{t,p,), xiS-iP) andu{t,p^) complex functions in real vari- 
able t and bar-homomorphic in complex variable p, (cf. (3.72)). Moreover, the above 
solution is smooth (analytic) if all a, p, k, cu, x o.i^'d ^ are smooth (analytic) functions. 
In particular, any solution of the equation (1.2) polynomial in x, y, z must be of the form 
(3.87) in which cr, p, uj, x o-nd v are polynomial in p,. 

Remcirk 3.4. In addition to the solutions in Theorems 3.1 and 3.2, the equation (1.2) 
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has the following simple solution: 

u = p,{V3{ay + bz)) x^, (3.88) 

where pi,{w) is the Weierstrass's elliptic function satisfying (2.39) and a, b are real con- 
stants such that + b"^ = 1. 

References 

[I] N. H. Ibragimov, Lie Group Analysis of Differential Equations, Volume 1, CRC Hand- 
book, CRC Press, 1995. 

[LRT] C. C. Lin, E. Reissner and H. S. Tsien, On two-dimensional non-steady motion 
of a slender body in a compressible fluid, J. Math, Phys. 27 (1948), no. 3, 220 

[Kpl] P. Kucharczyk, Group properties of the "short waves" equations in gas dynamics. 
Bull. Acad. Polon. Sci., Ser. Sci. Techn. XIII (1965), no. 5, 469 

[Kp2] P. Kucharczyk, Teoria Grup Liego w Zastosowaniu do Rowman Rozniczkowych 
Czaskowych, IPPT Polish Academy of Sciences, Warsaw, 65, 1967. 

[Ml] E. V. Mamontov, On the theory of nonstationary transonic flows, Dokl. Acad. 
Nauk SSSR 185 (1969), no. 3, 538 

[M2] E. V. Mamontov, Analytic perturbations in a nonstationary transonic stream, Di- 
namika Splosn. Sredy Vyp. 10 (1972), 217-222. 

[RS] O. S. Ryzhov and G. M. Shefter, On unsteady gas flows in Laval nozzles, Soviet 
Physics. Dokl. 4 (1959), 939-942. 

[Sg] G. D. Sevost'janov, An equation for nonstationary transonic flows of an ideal gas, 
Izv. Acad. Nauk SSSR Meh. Zidk. Gaza, 1977, no.l, 105-109. 

[Sv] S. V. Sukhinin, Group property and conservation laws of the equation of transonic 
motion of gas, Dinamika Splosh. Sredi 36 (1978), 130. 



19 



